This is an open assessment looking at potential health effects of a national fish promotion program in Finland. The details of the assessment are described at Opasnet. This file contains the R code to run the assessment model.
Knit to html for best performance.
Calculation is based on BONUS GOHERR project and its http://en.opasnet.org/w/Goherr_assessment.
What needs to be done for PFAS assessment?
objects.latest("Op_en2261", code_name="RR2") # RR on page [[Health impact assessment]]
## Loading objects:
## RR
RRorig <- RR@formula
RR@formula <- function(...) {
out <- RRorig()
out <- out * Ovariable(output=data.frame(Age=c(
"0 - 4", "5 - 9", "10 - 14", "15 - 19", "20 - 24", "25 - 29", "30 - 34", "35 - 39", "40 - 44", "45 - 49", "50 - 54",
"55 - 59", "60 - 64", "65 - 69", "70 - 74", "75 - 79", "80 - 84", "85+", "Undefined"),
Result=1),
marginal=c(TRUE,FALSE)
)
return(out)
}
# This code was forked from https://github.com/jtuomist/fishhealth/blob/master/fishhealth.Rmd
# This code was previously forked from code Op_fi5923/model on page [[Kotimaisen kalan edistämisohjelma]]
# The code was even more previously forked from Op_fi5889/model on page [[Ruori]] and Op_en7748/model on page [[Goherr assessment]]
dat <- opbase.data("Op_fi5932", subset="Malliparametrit")[-1] # [[PFAS-yhdeisteiden tautitaakka]]
dec <- opbase.data("Op_fi5932", subset="Decisions")[-1]
DecisionTableParser(dec)
CTable <- opbase.data("Op_fi5932",subset="CollapseMarginals")
#for(i in 1:ncol(CTable)) {CTable[[i]] <- as.character(CTable[[i]])} # The default is currently character, not factor
CollapseTableParser(CTable)
cat("Laskennassa käytetty data.\n")
## Laskennassa käytetty data.
dat
cat("Tarkastellut päätökset.\n")
## Tarkastellut päätökset.
dec
cat("Aggregoidut marginaalit.\n")
## Aggregoidut marginaalit.
CTable
dummy <- Ovariable("dummy",data=data.frame(Age="dummy", Fish="dummy", Compound="dummy", Area="dummy", Result=1)) # Keep these columns marginals
fish_proportion <- Ovariable( # How population subgroups eat fish differently
"fish_proportion",
data = data.frame(prepare(dat,"fish proportion",c("Type","Exposure_agent","Response","Unit"))),
unit="proportion of the mean")
total_amount <- Ovariable(
"total_amount",
data=prepare(dat, "amount", c("Type","Response","Exposure_agent","Unit")),
unit="M kg/a")
amount <- Ovariable(
"amount",
dependencies = data.frame(Name=c("total_amount", "fish_proportion")),
formula = function(...) {
amount <- total_amount
# Filleted weight, i.e. no loss.
amount <- amount * 1000 / 5.52 /365.25
# M kg/a per 5.52M population --> g/d per average person.
amount <- amount * fish_proportion
# fish_proportion tells the relative amount in each subgroup
# Match KKE-classification in amount with Fineli classification
tmp <- Ovariable(
output = data.frame(
Kala = c("Kasvatettu", "Kaupallinen", "Kirjolohi", "Silakka", "Vapaa-ajan", "Muu tuonti", "Tuontikirjolohi", "Tuontilohi"),
Fish = c("Whitefish", "Average fish","Rainbow trout", "Herring", "Average fish", "Average fish", "Rainbow trout", "Salmon"),
Result = 1
),
marginal = c(TRUE, TRUE, FALSE)
)
amount <- amount * tmp
return(amount)
},
unit="g/d"
)
# Exposure:To child and To eater not needed, because dioxins are not (yet) included
population <- Ovariable(
"population",
data=prepare(dat, "population", c("Type", "Exposure_agent", "Response","Unit")),
unit = "#"
)
population@data$Age <- factor(population@data$Age, levels=c(
"0 - 4", "5 - 9", "10 - 14", "15 - 19", "20 - 24", "25 - 29", "30 - 34", "35 - 39", "40 - 44", "45 - 49", "50 - 54",
"55 - 59", "60 - 64", "65 - 69", "70 - 74", "75 - 79", "80 - 84", "85+", "Undefined"))
incidence <- Ovariable(
"incidence",
data = prepare(dat,"incidence",c("Type","Exposure_agent","Unit")),
unit="1/person-year")
#incidence@data$Age[is.na(incidence@data$Age)] <- ""
case_burden <- Ovariable(
"case_burden",
data = prepare(dat,"case burden",c("Type", "Exposure_agent","Unit")),
unit="DALY/case")
ERFchoice <- Ovariable(
"ERFchoice",
data =
prepare(dat, "ERFchoice", c("Unit", "Type"))
)
InpBoD <- EvalOutput(Ovariable( # Evaluated because is not a dependency but an Input
"InpBoD",
data = prepare(dat, "BoD", c("Type","Exposure_agent","Unit")),
unit="DALY/a"
))
InpBoD$Response[InpBoD$Response=="All causes"] <- "All-cause mortality"
InpBoD$Response[InpBoD$Response=="Depressive disorders"] <- "Depression"
InpBoD$Response[InpBoD$Response=="Neoplasms"] <- "Cancer morbidity"
InpBoD$Response[InpBoD$Response=="Respiratory infections and tuberculosis"] <- "Immunosuppression" # Infections of 0-9-year-olds are assumed to represent the background BoD of immunosuppressive effect of PFAS
InpBoD$Response[InpBoD$Response=="Cardiovascular diseases"] <- "CHD2 mortality"
conc_vit <- Ovariable(
"conc_vit",
ddata = "Op_en1838", # [[Concentrations of beneficial nutrients in fish]]
subset = "Fineli data for common fish species",
unit="ALA mg/g, DHA mg/g, Fish g/g, Omega3 mg/g, Vitamin D µg/g f.w. after adjustment"
)
df = conc_vit@data
df$Nutrient[df$Nutrient=="D-vitamiini (µg)"] <- "Vitamin D"
df$Nutrient[df$Nutrient=="rasvahapot n-3 moni-tyydyttymättömät (g)"] <- "Omega3"
df$Nutrient[df$Nutrient=="rasvahappo 18:3 n-3 (alfalinoleenihappo) (mg)"] <- "ALA"
df$Nutrient[df$Nutrient=="rasvahappo 22:6 n-3 (DHA) (mg)"] <- "DHA"
df$Nutrient[df$Nutrient=="proteiini (g)"] <- "Fish"
df$conc_vitResult[df$Nutrient=="Fish"] <- "1"
df <- dropall(df[df$Nutrient %in% c("Vitamin D", "Omega3", "ALA", "DHA", "Fish") , ])
conc_vit@data <- df
######## Concentration of PFAS
# Data from EU-kalat3 (Finland excl Vanhankaupunginlahti): # pg/g fresh weight
# POP mean sd min Q0.025 median Q0.975 max
# 2.5% PFOS 2055.757 1404.045 305.0399 330.1365 1533.269 5029.697 5814.935
# Data from EU-kalat3 (Vanhankaupunginlahti, Helsinki) # ng/g f.w.
# POP mean sd min Q0.025 median Q0.975 max
#2.5% PFOS 14.428 11.94542 1.499441 1.607789 15.64988 35.32517 38.91994
conc_eukalat <- EvalOutput(Ovariable(
"conc_eukalat",
data = data.frame(
Area = c("Suomi","Helsinki"),
Compound="PFOS",
Result=c("2.056 (3.301 - 5.030)", "14.428 (1.499 - 35.325)")),
unit="ng/g fresh weight"
))
conc_pfas_raw <- Ovariable(
"conc_pfas_raw",
ddata="Op_fi5932",subset="PFAS concentrations",
unit="ng/g f.w.")
conc_pfas <- Ovariable(
"conc_pfas",
dependencies = data.frame(Name="conc_pfas_raw"),
formula = function(...) {
out <- conc_pfas_raw
out$Source <- NULL
out$Area <- "Porvoo"
out <- expand_index(out, list(Fish=list(Perch=c(
"Average fish", "Pike","Rainbow trout","Roach", "Salmon", "Vendace", "Whitefish"))))
out@marginal <- c(TRUE, TRUE, FALSE, TRUE, TRUE)
return(out)
},
unit="ng/g fresh weight"
)
#sum_pfas <- oapply(conc_pfas, cols=c("Kala","Compound"), FUN=sum)
#tmp <- conc_pfas / sum_pfas
#summary(tmp, marginals="Compound")
#
## This tells that PFOS consists of 71 - 97 % of the four key PFAS, while PFOA, PFNA, and PFHxS consist of
# 0 - 10 %, 2 - 18 %, and 0 - 9 %, respectively.
# Even if we included the next most abundant congeners, i.e. PFDA and PFUnA, the overall picture would not change.
conc <- Ovariable(
"conc",
dependencies = data.frame(
Name=c("conc_vit", "conc_pfas","conc_pcddf","conc_mehg"),
Ident = c(NA,NA,"Op_en3104/pcddf_mean","Op_en4004/conc_mehg_allfish")),
formula = function(...){
conc_vit <- oapply(conc_vit, cols=c("Kala", "Adjust"),FUN=mean)
colnames(conc_vit@output)[colnames(conc_vit@output)=="Nutrient"] <- "Compound"
conc_pfas <- oapply(conc_pfas, cols=c("Obs","Area"), FUN=mean)
conc_pfas$Compound[conc_pfas$Compound %in% c("PFOA","PFNA","PFHxS","PFOS")] <- "PFAS"
conc_pfas <- oapply(conc_pfas, cols="", FUN=sum)
conc_pcddf$Compound <- "TEQ"
conc_mehg <- oapply(conc_mehg[conc_mehg$Area=="Porvoo",], NULL, mean, c("Area","Kala")) # UPDATE FOR PROBABILISTIC APPROACH
conc_mehg$Compound <- "MeHg"
out <- OpasnetUtils::combine(conc_vit, conc_pfas, conc_pcddf, conc_mehg)
return(out)
}
)
#conc <- EvalOutput(conc, verbose=TRUE)
#View(conc@output)
#View(conc_pcddf@output)
#View(conc_mehg@output)
###################################################################
# Code copied from http://en.opasnet.org/w/Goherr_assessment#
mc2dparam<- list(
N2 = 10, # Number of iterations in the new Iter
strength = 50, # Sample size to which the fun is to be applied. Resembles number of observations
run2d = FALSE, # Should the mc2d function be used or not?
info = 1, # Ovariable that contains additional indices, e.g. newmarginals.
newmarginals = c("Exposure"), # Names of columns that are non-marginals but should be sampled enough to become marginals
method = "bootstrap", # which method to use for 2D Monte Carlo? Currently bootsrap is the only option.
fun = mean # Function for aggregating the first Iter dimension.
)
if(FALSE) {
## Exposure with background exposure but without mother's exposure to child
expo_dir <- Ovariable(
"expo_dir",
dependencies=data.frame(Name=c("amount","conc","expo_bg")),
formula = function(...) {
out <- conc[conc$Exposure_agent=="TEQ",] * 0 + 1
out$Exposure_agent <- "Fish"
out <- combine(conc, out, name="conc")
out <- oapply(amount * out, cols="Fish", FUN=sum)
out <- Ovariable(output = data.frame(
Exposcen = c("BAU", "No exposure"),
Result = c(1, 0)
), marginal=c(TRUE,FALSE)) * out + expo_bg
out$Exposure <- as.factor(
ifelse(
out$Exposure_agent %in% c("DHA", "MeHg"),
"To child",
"To eater"
)
)
return(out)
},
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d; Fish: g /d"
)
## Background-exposure to vitamin D and omega-3
addexposure <- Ovariable(
"addexposure",
ddata = "Op_en7748", # [[Benefit-risk assessment of Baltic herring and salmon intake]]
subset = "Background exposure",
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d"
)
# Should the background be specific for gender and country? At the moment it is.
expo_bg <- Ovariable(
"expo_bg",
dependencies = data.frame(Name="addexposure","info"),
formula = function(...) {
out <- addexposure
# Empty values ("") in indices must be replaced by NA so that Ops works correctly.
levels(out$Gender)[levels(out$Gender) == ""] <- NA
levels(out$Country)[levels(out$Country) == ""] <- NA
levels(out$Exposure_agent)[levels(out$Exposure_agent) == ""] <- NA
out@output <- fillna(out@output, c("Country", "Gender", "Exposure_agent"))
temp1 <- out[out$Exposure_agent %in% c("PCDDF","PCB") , ]
temp1 <- oapply(temp1, cols = "Exposure_agent", FUN = sum)
temp1$Exposure_agent <- "TEQ"
temp2 <- out[out$Exposure_agent %in% c("EPA", "DHA") , ]
temp2 <- oapply(temp2, cols = "Exposure_agent", FUN = sum)
temp2$Exposure_agent <- "Omega3"
out <- OpasnetUtils::combine(out, temp1, temp2)
out <- unkeep(out * info, prevresults = TRUE, sources = TRUE)
return(out)
},
unit = "PCDDF, PCB, TEQ: pg /d; Vitamin D, MeHg: µg /d; DHA, EPA, Omega3: mg /d"
)
# Stores non-marginal columns for further use.
info <- Ovariable(
"info",
dependencies = data.frame(Name = c("jsp")),
formula = function(...) {
out <- jsp
out$Group <- factor(
paste(out$Gender, out$Ages),
levels = c("Female 18-45", "Male 18-45", "Female >45", "Male >45")
)
out$Country <- factor(out$Country, ordered=FALSE)
out <- unique(out@output[c("Iter","Country","Group","Gender","Row")])
out$Result <- 1
return(out)
}
)
} # END IF
############################### Code from Goherr assessment ends
expo_dir <- Ovariable(
"expo_dir",
dependencies = data.frame(Name=c("conc", "amount")),
formula = function(...) {
conc$Fish[conc$Fish %in% c("Perch","Pike-perch","Eel","")] <- "Average fish"
conc <- oapply(conc, cols="", FUN=mean)
out <- conc * amount
# colnames(out@output)[colnames(out@output)=="Nutrient"] <- "Exposure_agent"
return(out)
},
unit="ng/d"
)
exposure <- Ovariable(
"exposure",
dependencies = data.frame(
Name = c(
"expo_dir", # direct exposure, i.e. the person eats or breaths the exposure agent themself
"expo_indir", # indirect exposure, i.e. the person (typically fetus or infant) is exposed via someone else (mother)
"mc2d" # 2D Monte Carlo function
),
Ident = c(
NA,
"Op_en7797/expo_indir2", # [[Infant's dioxin exposure]] # expo_indir
"Op_en7805/mc2d") # [[Two-dimensional Monte Carlo]]
),
formula = function(...) {
out <- OpasnetUtils::combine(expo_dir, expo_indir)
out <- unkeep(out, "Source.1", sources=TRUE)
out <- mc2d(out)
out$Exposure[is.na(out$Exposure)] <- "Direct"
return(out)
},
unit = "PCDDF, PCB, TEQ: (To eater: pg /day; to child: pg /g fat); Vitamin D, MeHg: µg /day; DHA, EPA, Omega3: mg /day"
)
exposure <- Ovariable(
"exposure",
dependencies=data.frame(Name=c("amount","conc")),
formula = function(...) {
out <- amount * conc
colnames(out@output)[colnames(out@output)=="Compound"] <- "Exposure_agent"
return(out)
}
)
objects.latest("Op_en2261",code_name="BoDattr2") # [[Health impact assessment]]
## Loading objects:
## BoDattr
tryCatch(BoDattr <- EvalOutput(BoDattr, verbose=TRUE))
## Evaluating BoDattr ...
##
## - BoD fetched successfully!
##
## - PAF fetched successfully!
## - Evaluating BoD ...
## - - Evaluating incidence ...
##
## done(0.01 secs)!
## - - Checking incidence marginals ... Response, Age, incidenceSource recognized as marginal(s).
## Loading required package: reshape2
##
## Attaching package: 'reshape2'
## The following object is masked from 'package:tidyr':
##
## smiths
## - - Processing incidence decisions ... done!
## - - Evaluating case_burden ...
##
## done(0 secs)!
## - - Checking case_burden marginals ... Response, case_burdenSource recognized as marginal(s).
## - - Processing case_burden marginal collapses ... done!
## - - Evaluating population ...
##
## done(0 secs)!
## - - Checking population marginals ... Gender, Age, populationSource recognized as marginal(s).
##
## - done(0.3 secs)!
## - Checking BoD marginals ... Response, Age, incidenceSource, Adjust, Gender, populationSource, BoDSource recognized as marginal(s).
## - Processing BoD inputs ... done!
## - Processing BoD marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying BoD, found NAs in indices: Adjust, InpBoDSource. They were
## automatically filled using fillna, which may result in a multiplied population.
## Please check your ovariable before using oapply.
## done!
## - Evaluating PAF ...
##
## - - dose fetched successfully!
##
## - - ERF fetched successfully!
##
## - - frexposed fetched successfully!
##
## - - P_illness fetched successfully!
##
## - - sumExposcen fetched successfully!
##
## - - mc2d fetched successfully!
## - - Evaluating dose ...
##
## - - - BW fetched successfully!
## - - - Evaluating exposure ...
## - - - - Evaluating amount ...
## - - - - - Evaluating total_amount ...
##
## done(0 secs)!
## - - - - - Checking total_amount marginals ... Kala, Scenario, total_amountSource recognized as marginal(s).
## - - - - - Evaluating fish_proportion ...
##
## done(0 secs)!
## - - - - - Checking fish_proportion marginals ... Gender, fish_proportionSource recognized as marginal(s).
##
## ---- done(0.12 secs)!
## - - - - Checking amount marginals ... Kala, Scenario, total_amountSource, Gender, fish_proportionSource, Fish, amountSource recognized as marginal(s).
## - - - - Evaluating conc ...
##
## - - - - - conc_pcddf fetched successfully!
##
## - - - - - conc_mehg fetched successfully!
## - - - - - Evaluating conc_vit ...
##
## done(0.02 secs)!
## - - - - - Checking conc_vit marginals ... Kala, Fish, Nutrient, conc_vitSource recognized as marginal(s).
## - - - - - Processing conc_vit decisions ... done!
## - - - - - Evaluating conc_pfas ...
## - - - - - - Evaluating conc_pfas_raw ...
##
## done(0.02 secs)!
## - - - - - - Checking conc_pfas_raw marginals ... Fish, Compound, conc_pfas_rawSource recognized as marginal(s).
##
## ----- done(0.09 secs)!
## - - - - - Checking conc_pfas marginals ... Fish, Compound, conc_pfas_rawSource, Area, conc_pfasSource recognized as marginal(s).
## - - - - - Evaluating conc_pcddf ...
##
## - - - - - - TEF fetched successfully!
## - - - - - - Evaluating TEF ...
##
## - - - - - - - TEFversion fetched successfully!
## - - - - - - - Evaluating TEFraw ...
##
## done(0.01 secs)!
## - - - - - - - Checking TEFraw marginals ... Group, Compound, TEFversion, TEFrawSource recognized as marginal(s).
##
## ------ done(27.63 secs)!
## - - - - - - Checking TEF marginals ... Compound, Group, TEFversion, TEFrawSource, TEFSource recognized as marginal(s).
##
## ----- done(1.03 mins)!
## - - - - - Checking conc_pcddf marginals ... Fish, conc_pcddfSource recognized as marginal(s).
## - - - - - Evaluating conc_mehg ...
##
## - - - - - - Hg fetched successfully!
##
## ----- done(23.03 secs)!
## - - - - - Checking conc_mehg marginals ... Area, Kala, Fish, conc_mehgSource recognized as marginal(s).
##
## ---- done(2.22 mins)!
## - - - - Checking conc marginals ... Fish, Compound, conc_pfas_rawSource, concSource recognized as marginal(s).
##
## --- done(2.22 mins)!
## - - - Checking exposure marginals ... Kala, Scenario, total_amountSource, Gender, fish_proportionSource, Fish, amountSource, Exposure_agent, conc_pfas_rawSource, concSource, exposureSource recognized as marginal(s).
## - - - Processing exposure marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying exposure, found NAs in indices: conc_pfas_rawSource. They were
## automatically filled using fillna, which may result in a multiplied population.
## Please check your ovariable before using oapply.
## done!
## - - - Evaluating BW ...
##
## done(0 secs)!
## - - - Checking BW marginals ... BWSource recognized as marginal(s).
##
## -- done(2.69 mins)!
## - - Checking dose marginals ... Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, Exposure_agent, conc_pfas_rawSource, concSource, Scaling, exposureSource, BWSource, doseSource recognized as marginal(s).
## - - Processing dose marginal collapses ... done!
## - - Evaluating ERF ...
##
## - - - ERF_env fetched successfully!
##
## - - - ERF_omega3 fetched successfully!
##
## - - - ERF_mehg fetched successfully!
##
## - - - ERF_diox fetched successfully!
##
## - - - ERF_vit fetched successfully!
##
## - - - ERF_micr fetched successfully!
##
## - - - ERF_pfas fetched successfully!
## - - - Evaluating ERF_env ...
##
## done(0.02 secs)!
## - - - Checking ERF_env marginals ... Exposure_agent, Response, Subgroup, Exposure, ER_function, Scaling, Exposure_unit, Observation, ERF_envSource recognized as marginal(s).
## - - - Evaluating ERF_omega3 ...
##
## done(0.01 secs)!
## - - - Checking ERF_omega3 marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_omega3Source recognized as marginal(s).
## - - - Evaluating ERF_mehg ...
##
## done(0 secs)!
## - - - Checking ERF_mehg marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_mehgSource recognized as marginal(s).
## - - - Evaluating ERF_diox ...
##
## done(0 secs)!
## - - - Checking ERF_diox marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_dioxSource recognized as marginal(s).
## - - - Evaluating ERF_vit ...
##
## done(0 secs)!
## - - - Checking ERF_vit marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_vitSource recognized as marginal(s).
## - - - Evaluating ERF_micr ...
##
## done(0 secs)!
## - - - Checking ERF_micr marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_micrSource recognized as marginal(s).
## - - - Evaluating ERF_pfas ...
##
## done(0 secs)!
## - - - Checking ERF_pfas marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_pfasSource recognized as marginal(s).
## - - - Evaluating ERFchoice ...
##
## done(0 secs)!
## - - - Checking ERFchoice marginals ... Exposure_agent, Response, Scaling, Exposure, ER_function, ERFchoiceSource recognized as marginal(s).
##
## -- done(2.95 mins)!
## - - Checking ERF marginals ... Exposure_agent, Response, Exposure, ER_function, Scaling, Observation, ERFSource recognized as marginal(s).
## - - Processing ERF marginal collapses ... done!
## - - Evaluating RR ...
## - - - Processing dose marginal collapses ... done!
## - - - Processing ERF marginal collapses ... done!
##
## -- done(0.18 secs)!
## - - Checking RR marginals ... Exposure_agent, Response, ER_function, Scaling, ERFSource, Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, conc_pfas_rawSource, doseSource, Age, RRSource recognized as marginal(s).
## - - Evaluating frexposed ...
##
## done(0 secs)!
## - - Checking frexposed marginals ... frexposedSource recognized as marginal(s).
## - - Evaluating P_illness ...
##
## done(0 secs)!
## - - Checking P_illness marginals ... Response, Illness, Age, P_illnessSource recognized as marginal(s).
##
## - done(8.1 mins)!
## - Checking PAF marginals ... Exposure_agent, Response, ER_function, Scaling, ERFSource, Scenario, total_amountSource, Gender, fish_proportionSource, amountSource, conc_pfas_rawSource, doseSource, frexposedSource, Age, incidenceSource, Adjust, RRSource, PAFSource recognized as marginal(s).
## - Processing PAF marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying PAF, found NAs in indices: Age, Adjust. They were automatically
## filled using fillna, which may result in a multiplied population. Please check
## your ovariable before using oapply.
## done!
##
## done(8.91 mins)!
## Checking BoDattr marginals ... Response, Age, Gender, Adjust, InpBoDSource, Exposure_agent, Scenario, total_amountSource, fish_proportionSource, amountSource, conc_pfas_rawSource, PAFSource, BoDattrSource recognized as marginal(s).
oprint(summary(amount,"mean"))
## Kala Scenario Gender Fish mean
## 1 Kaupallinen BAU Female Average fish 1.2697279
## 2 Muu tuonti BAU Female Average fish 11.1101191
## 3 Vapaa-ajan BAU Female Average fish 3.9282207
## 4 Kaupallinen BAU Male Average fish 1.9045919
## 5 Muu tuonti BAU Male Average fish 16.6651787
## 6 Vapaa-ajan BAU Male Average fish 5.8923310
## 7 Silakka BAU Female Herring 0.6348640
## 8 Silakka BAU Male Herring 0.9522959
## 9 Kirjolohi BAU Female Rainbow trout 2.6584928
## 10 Tuontikirjolohi BAU Female Rainbow trout 1.9045919
## 11 Kirjolohi BAU Male Rainbow trout 3.9877392
## 12 Tuontikirjolohi BAU Male Rainbow trout 2.8568878
## 13 Tuontilohi BAU Female Salmon 9.3642433
## 14 Tuontilohi BAU Male Salmon 14.0463649
## 15 Kasvatettu BAU Female Whitefish 0.2380740
## 16 Kasvatettu BAU Male Whitefish 0.3571110
oprint(summary(BoD,marginals=c("Age","Response"),"mean"))
## Age Response mean
## 1 0 - 4 Cancer morbidity 309.5200
## 2 10 - 14 Cancer morbidity 347.0750
## 3 15 - 19 Cancer morbidity 414.9750
## 4 20 - 24 Cancer morbidity 569.5250
## 5 25 - 29 Cancer morbidity 795.4550
## 6 30 - 34 Cancer morbidity 1147.9850
## 7 35 - 39 Cancer morbidity 1679.1250
## 8 40 - 44 Cancer morbidity 2452.5650
## 9 45 - 49 Cancer morbidity 3849.3800
## 10 5 - 9 Cancer morbidity 345.5200
## 11 50 - 54 Cancer morbidity 6987.9000
## 12 55 - 59 Cancer morbidity 11237.3300
## 13 60 - 64 Cancer morbidity 16331.9100
## 14 65 - 69 Cancer morbidity 21899.6500
## 15 70 - 74 Cancer morbidity 25551.8600
## 16 75 - 79 Cancer morbidity 18015.2300
## 17 80 - 84 Cancer morbidity 14650.9150
## 18 85 - 89 Cancer morbidity 8413.0550
## 19 90 - 94 Cancer morbidity 3537.1750
## 20 Undefined Dioxin recommendation tolerable daily intake 305.4067
## 21 Undefined Dioxin recommendation tolerable daily intake 2018 888.4560
## 22 0 - 4 Immunosuppression 295.9850
## 23 5 - 9 Immunosuppression 261.5350
## 24 0 - 4 Loss in child's IQ points 16778.3774
## 25 Undefined PFAS TWI 2776.4249
## 26 0 - 4 Sperm concentration 4478.6700
## 27 Undefined Vitamin D recommendation 610.8135
## 28 0 - 4 Yes or no dental defect 343.9619
## 29 0 - 4 All-cause mortality 4867.4650
## 30 10 - 14 All-cause mortality 1085.5450
## 31 15 - 19 All-cause mortality 2965.9450
## 32 20 - 24 All-cause mortality 5274.8450
## 33 25 - 29 All-cause mortality 6263.5400
## 34 30 - 34 All-cause mortality 6710.3100
## 35 35 - 39 All-cause mortality 8275.6450
## 36 40 - 44 All-cause mortality 10232.0150
## 37 45 - 49 All-cause mortality 13691.2800
## 38 5 - 9 All-cause mortality 931.1250
## 39 50 - 54 All-cause mortality 21321.8650
## 40 55 - 59 All-cause mortality 29861.9100
## 41 60 - 64 All-cause mortality 39200.8750
## 42 65 - 69 All-cause mortality 49295.5050
## 43 70 - 74 All-cause mortality 59799.7000
## 44 75 - 79 All-cause mortality 49428.5150
## 45 80 - 84 All-cause mortality 51968.9150
## 46 85 - 89 All-cause mortality 42765.3450
## 47 90 - 94 All-cause mortality 26756.6400
## 48 0 - 4 Depression 0.4150
## 49 10 - 14 Depression 555.8850
## 50 15 - 19 Depression 1320.8250
## 51 20 - 24 Depression 1841.0150
## 52 25 - 29 Depression 1801.4650
## 53 30 - 34 Depression 1590.0200
## 54 35 - 39 Depression 1721.7250
## 55 40 - 44 Depression 1700.2800
## 56 45 - 49 Depression 1534.1700
## 57 5 - 9 Depression 76.4850
## 58 50 - 54 Depression 1718.9550
## 59 55 - 59 Depression 1767.1650
## 60 60 - 64 Depression 1676.7400
## 61 65 - 69 Depression 1561.5750
## 62 70 - 74 Depression 1426.2350
## 63 75 - 79 Depression 799.7000
## 64 80 - 84 Depression 581.2950
## 65 85 - 89 Depression 363.0750
## 66 90 - 94 Depression 186.1150
## 67 0 - 4 CHD2 mortality 53.6100
## 68 10 - 14 CHD2 mortality 75.3150
## 69 15 - 19 CHD2 mortality 145.3450
## 70 20 - 24 CHD2 mortality 254.1000
## 71 25 - 29 CHD2 mortality 437.6450
## 72 30 - 34 CHD2 mortality 636.4200
## 73 35 - 39 CHD2 mortality 1168.2600
## 74 40 - 44 CHD2 mortality 1977.5650
## 75 45 - 49 CHD2 mortality 3117.3100
## 76 5 - 9 CHD2 mortality 48.7800
## 77 50 - 54 CHD2 mortality 5476.2900
## 78 55 - 59 CHD2 mortality 8747.4650
## 79 60 - 64 CHD2 mortality 13077.9700
## 80 65 - 69 CHD2 mortality 18100.7500
## 81 70 - 74 CHD2 mortality 24549.2850
## 82 75 - 79 CHD2 mortality 23413.4850
## 83 80 - 84 CHD2 mortality 27899.3350
## 84 85 - 89 CHD2 mortality 25240.0350
## 85 90 - 94 CHD2 mortality 16630.6150
## 86 15 - 19 Breast cancer 2.0900
## 87 20 - 24 Breast cancer 7.7850
## 88 25 - 29 Breast cancer 36.8200
## 89 30 - 34 Breast cancer 145.1600
## 90 35 - 39 Breast cancer 263.6000
## 91 40 - 44 Breast cancer 433.3500
## 92 45 - 49 Breast cancer 704.9550
## 93 50 - 54 Breast cancer 1045.3250
## 94 55 - 59 Breast cancer 1341.8150
## 95 60 - 64 Breast cancer 1513.2150
## 96 65 - 69 Breast cancer 1649.3300
## 97 70 - 74 Breast cancer 1678.0350
## 98 75 - 79 Breast cancer 1232.0100
## 99 80 - 84 Breast cancer 950.8800
## 100 85 - 89 Breast cancer 578.1900
## 101 90 - 94 Breast cancer 296.1800
oprint(summary(BoDattr,marginals=c("Age","Response"),"mean"))
## Age Response mean
## 1 0 - 4 Cancer morbidity 8.092797e+01
## 2 10 - 14 Cancer morbidity 8.850176e+01
## 3 15 - 19 Cancer morbidity 1.079079e+02
## 4 20 - 24 Cancer morbidity 1.483914e+02
## 5 25 - 29 Cancer morbidity 2.074377e+02
## 6 30 - 34 Cancer morbidity 2.931146e+02
## 7 35 - 39 Cancer morbidity 4.269816e+02
## 8 40 - 44 Cancer morbidity 6.183735e+02
## 9 45 - 49 Cancer morbidity 9.690704e+02
## 10 5 - 9 Cancer morbidity 8.825348e+01
## 11 50 - 54 Cancer morbidity 1.781467e+03
## 12 55 - 59 Cancer morbidity 2.926500e+03
## 13 60 - 64 Cancer morbidity 4.292083e+03
## 14 65 - 69 Cancer morbidity 5.788529e+03
## 15 70 - 74 Cancer morbidity 6.750697e+03
## 16 75 - 79 Cancer morbidity 4.726815e+03
## 17 80 - 84 Cancer morbidity 3.804571e+03
## 18 Undefined Dioxin recommendation tolerable daily intake 3.054067e+02
## 19 Undefined Dioxin recommendation tolerable daily intake 2018 8.884560e+02
## 20 0 - 4 Immunosuppression 6.818641e+00
## 21 5 - 9 Immunosuppression 7.294478e+00
## 22 0 - 4 Loss in child's IQ points 6.974660e+02
## 23 Undefined PFAS TWI 1.370935e+03
## 24 0 - 4 Sperm concentration 4.548773e+03
## 25 Undefined Vitamin D recommendation 6.108135e+02
## 26 0 - 4 Yes or no dental defect 2.530885e+03
## 27 0 - 4 All-cause mortality -3.940215e+02
## 28 10 - 14 All-cause mortality -8.729181e+01
## 29 15 - 19 All-cause mortality -2.510075e+02
## 30 20 - 24 All-cause mortality -4.582433e+02
## 31 25 - 29 All-cause mortality -5.472993e+02
## 32 30 - 34 All-cause mortality -5.805901e+02
## 33 35 - 39 All-cause mortality -7.164739e+02
## 34 40 - 44 All-cause mortality -8.793680e+02
## 35 45 - 49 All-cause mortality -1.166310e+03
## 36 5 - 9 All-cause mortality -7.424681e+01
## 37 50 - 54 All-cause mortality -1.808052e+03
## 38 55 - 59 All-cause mortality -2.535812e+03
## 39 60 - 64 All-cause mortality -3.319267e+03
## 40 65 - 69 All-cause mortality -4.145066e+03
## 41 70 - 74 All-cause mortality -4.980650e+03
## 42 75 - 79 All-cause mortality -4.044850e+03
## 43 80 - 84 All-cause mortality -4.134937e+03
## 44 0 - 4 Depression -7.647127e-02
## 45 10 - 14 Depression -9.937062e+01
## 46 15 - 19 Depression -2.347323e+02
## 47 20 - 24 Depression -3.276647e+02
## 48 25 - 29 Depression -3.226701e+02
## 49 30 - 34 Depression -2.865779e+02
## 50 35 - 39 Depression -3.110681e+02
## 51 40 - 44 Depression -3.073361e+02
## 52 45 - 49 Depression -2.766914e+02
## 53 5 - 9 Depression -1.387565e+01
## 54 50 - 54 Depression -3.090770e+02
## 55 55 - 59 Depression -3.166728e+02
## 56 60 - 64 Depression -2.997941e+02
## 57 65 - 69 Depression -2.791761e+02
## 58 70 - 74 Depression -2.551780e+02
## 59 75 - 79 Depression -1.425020e+02
## 60 80 - 84 Depression -1.024717e+02
## 61 0 - 4 CHD2 mortality -8.359787e+00
## 62 10 - 14 CHD2 mortality -1.170566e+01
## 63 15 - 19 CHD2 mortality -2.262663e+01
## 64 20 - 24 CHD2 mortality -3.966234e+01
## 65 25 - 29 CHD2 mortality -6.844631e+01
## 66 30 - 34 CHD2 mortality -9.964232e+01
## 67 35 - 39 CHD2 mortality -1.833859e+02
## 68 40 - 44 CHD2 mortality -3.106654e+02
## 69 45 - 49 CHD2 mortality -4.899542e+02
## 70 5 - 9 CHD2 mortality -7.579246e+00
## 71 50 - 54 CHD2 mortality -8.613391e+02
## 72 55 - 59 CHD2 mortality -1.375930e+03
## 73 60 - 64 CHD2 mortality -2.055932e+03
## 74 65 - 69 CHD2 mortality -2.841726e+03
## 75 70 - 74 CHD2 mortality -3.847268e+03
## 76 75 - 79 CHD2 mortality -3.660036e+03
## 77 80 - 84 CHD2 mortality -4.346514e+03
## 78 15 - 19 Breast cancer -4.624153e-01
## 79 20 - 24 Breast cancer -1.587543e+00
## 80 25 - 29 Breast cancer -7.326368e+00
## 81 30 - 34 Breast cancer -2.873340e+01
## 82 35 - 39 Breast cancer -5.218217e+01
## 83 40 - 44 Breast cancer -8.577757e+01
## 84 45 - 49 Breast cancer -1.395958e+02
## 85 50 - 54 Breast cancer -2.070689e+02
## 86 55 - 59 Breast cancer -2.658472e+02
## 87 60 - 64 Breast cancer -2.998331e+02
## 88 65 - 69 Breast cancer -3.267095e+02
## 89 70 - 74 Breast cancer -3.327787e+02
## 90 75 - 79 Breast cancer -2.441402e+02
## 91 80 - 84 Breast cancer -1.884558e+02
oprint(summary(BoDattr,marginals=c("Exposure_agent","Response"),"mean"))
## Exposure_agent Response
## 1 TEQ Cancer morbidity
## 2 TEQ Dioxin recommendation tolerable daily intake
## 3 TEQ Dioxin recommendation tolerable daily intake 2018
## 4 PFAS Immunosuppression
## 5 DHA Loss in child's IQ points
## 6 MeHg Loss in child's IQ points
## 7 PFAS PFAS TWI
## 8 TEQ Sperm concentration
## 9 Vitamin D Vitamin D recommendation
## 10 TEQ Yes or no dental defect
## 11 Fish All-cause mortality
## 12 Fish Depression
## 13 Omega3 CHD2 mortality
## 14 Omega3 Breast cancer
## mean
## 1 1947.036623
## 2 305.406736
## 3 888.455958
## 4 7.056559
## 5 -3292.013959
## 6 4686.945878
## 7 1370.934688
## 8 4548.773342
## 9 610.813471
## 10 2530.884510
## 11 -1771.969747
## 12 -228.525606
## 13 -1190.045495
## 14 -155.749901
oprint(summary(BoDattr,marginals=c("Gender","Response"),"mean"))
## Gender Response mean
## 1 Female Cancer morbidity 1378.857757
## 2 Male Cancer morbidity 2515.215489
## 3 Female Dioxin recommendation tolerable daily intake 309.207840
## 4 Male Dioxin recommendation tolerable daily intake 301.605631
## 5 Female Dioxin recommendation tolerable daily intake 2018 899.513715
## 6 Male Dioxin recommendation tolerable daily intake 2018 877.398200
## 7 Female Immunosuppression 5.205862
## 8 Male Immunosuppression 8.907256
## 9 Female Loss in child's IQ points 542.752772
## 10 Male Loss in child's IQ points 852.179148
## 11 Female PFAS TWI 0.000000
## 12 Male PFAS TWI 2741.869376
## 13 Female Sperm concentration 3539.756033
## 14 Male Sperm concentration 5557.790651
## 15 Female Vitamin D recommendation 618.415679
## 16 Male Vitamin D recommendation 603.211263
## 17 Female Yes or no dental defect 1969.479031
## 18 Male Yes or no dental defect 3092.289990
## 19 Female All-cause mortality -1003.248948
## 20 Male All-cause mortality -2540.690546
## 21 Female Depression -235.343301
## 22 Male Depression -221.707910
## 23 Female CHD2 mortality -815.643510
## 24 Male CHD2 mortality -1564.447480
## 25 Female Breast cancer -308.799252
## 26 Male Breast cancer -2.700550
oprint(summary(case_burden,"mean"))
## Response mean
## 1 Dioxin recommendation tolerable daily intake 0.001004988
## 2 Dioxin recommendation tolerable daily intake 2018 0.001004988
## 3 Loss in child's IQ points 0.110000000
## 4 PFAS TWI 0.001004988
## 5 Sperm concentration 2.500000000
## 6 Vitamin D recommendation 0.001004988
## 7 Yes or no dental defect 0.060000000
oprint(summary(conc,"mean"))
## Fish Compound mean
## 1 Average fish ALA 0.69000000
## 2 Bream ALA 0.22000000
## 3 Herring ALA 1.74000000
## 4 Pike ALA 0.08000000
## 5 Rainbow trout ALA 4.81000000
## 6 Roach ALA 0.10000000
## 7 Salmon ALA 7.96000000
## 8 Vendace ALA 1.35000000
## 9 Whitefish ALA 2.22000000
## 10 Average fish DHA 2.54000000
## 11 Bream DHA 2.73000000
## 12 Herring DHA 5.86000000
## 13 Pike DHA 0.30000000
## 14 Rainbow trout DHA 7.57000000
## 15 Roach DHA 2.87000000
## 16 Salmon DHA 6.69000000
## 17 Vendace DHA 3.00000000
## 18 Whitefish DHA 3.94000000
## 19 Average fish Fish 1.00000000
## 20 Bream Fish 1.00000000
## 21 Herring Fish 1.00000000
## 22 Pike Fish 1.00000000
## 23 Rainbow trout Fish 1.00000000
## 24 Roach Fish 1.00000000
## 25 Salmon Fish 1.00000000
## 26 Vendace Fish 1.00000000
## 27 Whitefish Fish 1.00000000
## 28 Average fish MeHg 0.11413395
## 29 Bream MeHg 0.63339772
## 30 Burbot MeHg 1.02686230
## 31 Chaetiliidae MeHg 0.01634649
## 32 Flounder MeHg 0.04582576
## 33 Fourhorn sculpin MeHg 0.64990137
## 34 Herring MeHg 0.02680994
## 35 Perch MeHg 0.14874312
## 36 Pike MeHg 0.90310953
## 37 Pike-perch MeHg 0.25432475
## 38 Roach MeHg 0.26739658
## 39 Sprat MeHg 0.03187976
## 40 White bream MeHg 1.21073890
## 41 Whitefish MeHg 0.06619502
## 42 Average fish Omega3 7.00000000
## 43 Bream Omega3 6.00000000
## 44 Herring Omega3 24.00000000
## 45 Pike Omega3 0.50000000
## 46 Rainbow trout Omega3 18.00000000
## 47 Roach Omega3 5.00000000
## 48 Salmon Omega3 23.00000000
## 49 Vendace Omega3 10.00000000
## 50 Whitefish Omega3 10.00000000
## 51 Average fish PFAS 8.79312500
## 52 Bream PFAS 4.95000000
## 53 Eel PFAS 7.89500000
## 54 Herring PFAS 1.65500000
## 55 Perch PFAS 8.79312500
## 56 Pike PFAS 8.79312500
## 57 Pike-perch PFAS 2.59500000
## 58 Rainbow trout PFAS 8.79312500
## 59 Roach PFAS 8.79312500
## 60 Salmon PFAS 8.79312500
## 61 Vendace PFAS 8.79312500
## 62 Whitefish PFAS 8.79312500
## 63 Baltic herring TEQ 90.08722156
## 64 Bream TEQ 6.20175433
## 65 Burbot TEQ 0.82549312
## 66 Cod TEQ 1.67367981
## 67 Flounder TEQ 15.36558024
## 68 Perch TEQ 3.56897420
## 69 Pike TEQ 2.14862962
## 70 Pike-perch TEQ 2.70814568
## 71 River lamprey TEQ 7.34301416
## 72 Roach TEQ 0.72124877
## 73 Salmon TEQ 20.00778296
## 74 Sea trout TEQ 9.31778459
## 75 Sprat TEQ 4.15456559
## 76 Vendace TEQ 0.99253888
## 77 Whitefish TEQ 5.74967872
## 78 Average fish Vitamin D 0.10500000
## 79 Bream Vitamin D 0.14000000
## 80 Herring Vitamin D 0.15600000
## 81 Pike Vitamin D 0.02100000
## 82 Rainbow trout Vitamin D 0.05100000
## 83 Roach Vitamin D 0.10000000
## 84 Salmon Vitamin D 0.06700000
## 85 Vendace Vitamin D 0.09400000
## 86 Whitefish Vitamin D 0.14400000
oprint(summary(dose,"mean"))
## Scenario Gender Exposure_agent Scaling mean
## 1 BAU Female ALA BW 1.56247954
## 2 BAU Male ALA BW 2.34371931
## 3 BAU Female DHA BW 2.04671635
## 4 BAU Male DHA BW 3.07007453
## 5 BAU Female Fish BW 0.44440477
## 6 BAU Male Fish BW 0.66660715
## 7 BAU Female MeHg BW 0.02705835
## 8 BAU Male MeHg BW 0.04058752
## 9 BAU Female Omega3 BW 6.13267239
## 10 BAU Male Omega3 BW 9.19900859
## 11 BAU Female PFAS BW 3.84296754
## 12 BAU Male PFAS BW 5.76445130
## 13 BAU Female TEQ BW 2.69609423
## 14 BAU Male TEQ BW 4.04414134
## 15 BAU Female Vitamin D BW 0.03865415
## 16 BAU Male Vitamin D BW 0.05798122
## 17 BAU Female ALA Log10 2.03891238
## 18 BAU Male ALA Log10 2.21500364
## 19 BAU Female DHA Log10 2.15615570
## 20 BAU Male DHA Log10 2.33224696
## 21 BAU Female Fish Log10 1.49287675
## 22 BAU Male Fish Log10 1.66896801
## 23 BAU Female MeHg Log10 0.27739929
## 24 BAU Male MeHg Log10 0.45349055
## 25 BAU Female Omega3 Log10 2.63274781
## 26 BAU Male Omega3 Log10 2.80883906
## 27 BAU Female PFAS Log10 2.42976476
## 28 BAU Male PFAS Log10 2.60585601
## 29 BAU Female TEQ Log10 2.27583311
## 30 BAU Male TEQ Log10 2.45192437
## 31 BAU Female Vitamin D Log10 0.43229411
## 32 BAU Male Vitamin D Log10 0.60838537
## 33 BAU Female ALA None 109.37356784
## 34 BAU Male ALA None 164.06035175
## 35 BAU Female DHA None 143.27014453
## 36 BAU Male DHA None 214.90521680
## 37 BAU Female Fish None 31.10833358
## 38 BAU Male Fish None 46.66250037
## 39 BAU Female MeHg None 1.89408422
## 40 BAU Male MeHg None 2.84112633
## 41 BAU Female Omega3 None 429.28706762
## 42 BAU Male Omega3 None 643.93060143
## 43 BAU Female PFAS None 269.00772748
## 44 BAU Male PFAS None 403.51159123
## 45 BAU Female TEQ None 188.72659592
## 46 BAU Male TEQ None 283.08989388
## 47 BAU Female Vitamin D None 2.70579016
## 48 BAU Male Vitamin D None 4.05868524
oprint(summary(ERF,"mean"))
## Exposure_agent Response
## 1 Fish All-cause mortality
## 2 Omega3 Breast cancer
## 3 Fish Depression
## 4 DHA Loss in child's IQ points
## 5 TEQ Sperm concentration
## 6 TEQ Yes or no dental defect
## 7 Omega3 CHD2 mortality
## 8 Vitamin D Vitamin D recommendation
## 9 MeHg Loss in child's IQ points
## 10 PFAS Immunosuppression
## 11 PFAS PFAS TWI
## 12 TEQ Cancer morbidity
## 13 TEQ Dioxin recommendation tolerable daily intake
## 14 TEQ Dioxin recommendation tolerable daily intake 2018
## 15 Fish All-cause mortality
## 16 Omega3 Breast cancer
## 17 Fish Depression
## 18 DHA Loss in child's IQ points
## 19 TEQ Sperm concentration
## 20 TEQ Yes or no dental defect
## 21 Omega3 CHD2 mortality
## 22 Vitamin D Vitamin D recommendation
## 23 MeHg Loss in child's IQ points
## 24 PFAS Immunosuppression
## 25 PFAS PFAS TWI
## 26 TEQ Cancer morbidity
## 27 TEQ Dioxin recommendation tolerable daily intake
## 28 TEQ Dioxin recommendation tolerable daily intake 2018
## ER_function Scaling Observation mean
## 1 RR None ERF 0.997871700
## 2 RR None ERF 0.999487200
## 3 RR None ERF 0.994690400
## 4 ERS None ERF -0.001300000
## 5 ERS None ERF 0.000060000
## 6 ERS None ERF 0.001390971
## 7 Relative Hill None ERF -0.170000000
## 8 Step None ERF 100.000000000
## 9 ERS BW ERF 9.800000000
## 10 ERS BW ERF 0.022700000
## 11 TWI BW ERF 4.400000000
## 12 CSF BW ERF 0.000500000
## 13 TDI BW ERF 2.000000000
## 14 TDI BW ERF 0.288900000
## 15 RR None Threshold 0.000000000
## 16 RR None Threshold 0.000000000
## 17 RR None Threshold 0.000000000
## 18 ERS None Threshold 0.000000000
## 19 ERS None Threshold 0.000000000
## 20 ERS None Threshold 0.000000000
## 21 Relative Hill None Threshold 47.000000000
## 22 Step None Threshold 10.000000000
## 23 ERS BW Threshold 0.000000000
## 24 ERS BW Threshold 0.000000000
## 25 TWI BW Threshold 0.000000000
## 26 CSF BW Threshold 0.000000000
## 27 TDI BW Threshold 0.000000000
## 28 TDI BW Threshold 0.000000000
oprint(summary(expo_dir,"mean"))
## Warning in oapply(conc, cols = "", FUN = mean): While oapplying conc, found NAs
## in indices: conc_pfas_rawSource. They were automatically filled using fillna,
## which may result in a multiplied population. Please check your ovariable before
## using oapply.
## Fish Compound Kala Scenario Gender mean
## 1 Whitefish ALA Kasvatettu BAU Female 0.52852424
## 2 Whitefish DHA Kasvatettu BAU Female 0.93801149
## 3 Whitefish Fish Kasvatettu BAU Female 0.23807398
## 4 Whitefish MeHg Kasvatettu BAU Female 0.01575931
## 5 Whitefish Omega3 Kasvatettu BAU Female 2.38073981
## 6 Whitefish PFAS Kasvatettu BAU Female 2.09341428
## 7 Whitefish TEQ Kasvatettu BAU Female 1.36884891
## 8 Whitefish Vitamin D Kasvatettu BAU Female 0.03428265
## 9 Average fish ALA Kaupallinen BAU Female 0.87611225
## 10 Average fish DHA Kaupallinen BAU Female 3.22510887
## 11 Average fish Fish Kaupallinen BAU Female 1.26972790
## 12 Average fish MeHg Kaupallinen BAU Female 0.21890186
## 13 Average fish Omega3 Kaupallinen BAU Female 8.88809531
## 14 Average fish PFAS Kaupallinen BAU Female 8.91229950
## 15 Average fish TEQ Kaupallinen BAU Female 3.98511713
## 16 Average fish Vitamin D Kaupallinen BAU Female 0.13332143
## 17 Rainbow trout ALA Kirjolohi BAU Female 12.78735034
## 18 Rainbow trout DHA Kirjolohi BAU Female 20.12479045
## 19 Rainbow trout Fish Kirjolohi BAU Female 2.65849279
## 20 Rainbow trout Omega3 Kirjolohi BAU Female 47.85287028
## 21 Rainbow trout PFAS Kirjolohi BAU Female 23.37645944
## 22 Rainbow trout Vitamin D Kirjolohi BAU Female 0.13558313
## 23 Average fish ALA Muu tuonti BAU Female 7.66598220
## 24 Average fish DHA Muu tuonti BAU Female 28.21970261
## 25 Average fish Fish Muu tuonti BAU Female 11.11011914
## 26 Average fish MeHg Muu tuonti BAU Female 1.91539131
## 27 Average fish Omega3 Muu tuonti BAU Female 77.77083395
## 28 Average fish PFAS Muu tuonti BAU Female 77.98262060
## 29 Average fish TEQ Muu tuonti BAU Female 34.86977486
## 30 Average fish Vitamin D Muu tuonti BAU Female 1.16656251
## 31 Herring ALA Silakka BAU Female 1.10466327
## 32 Herring DHA Silakka BAU Female 3.72030275
## 33 Herring Fish Silakka BAU Female 0.63486395
## 34 Herring MeHg Silakka BAU Female 0.01702066
## 35 Herring Omega3 Silakka BAU Female 15.23673482
## 36 Herring PFAS Silakka BAU Female 1.05069984
## 37 Herring Vitamin D Silakka BAU Female 0.09903878
## 38 Rainbow trout ALA Tuontikirjolohi BAU Female 9.16108681
## 39 Rainbow trout DHA Tuontikirjolohi BAU Female 14.41776032
## 40 Rainbow trout Fish Tuontikirjolohi BAU Female 1.90459185
## 41 Rainbow trout Omega3 Tuontikirjolohi BAU Female 34.28265333
## 42 Rainbow trout PFAS Tuontikirjolohi BAU Female 16.74731423
## 43 Rainbow trout Vitamin D Tuontikirjolohi BAU Female 0.09713418
## 44 Salmon ALA Tuontilohi BAU Female 74.53937644
## 45 Salmon DHA Tuontilohi BAU Female 62.64678749
## 46 Salmon Fish Tuontilohi BAU Female 9.36424327
## 47 Salmon Omega3 Tuontilohi BAU Female 215.37759525
## 48 Salmon PFAS Tuontilohi BAU Female 82.34096162
## 49 Salmon TEQ Tuontilohi BAU Female 187.35774701
## 50 Salmon Vitamin D Tuontilohi BAU Female 0.62740430
## 51 Average fish ALA Vapaa-ajan BAU Female 2.71047228
## 52 Average fish DHA Vapaa-ajan BAU Female 9.97768056
## 53 Average fish Fish Vapaa-ajan BAU Female 3.92822069
## 54 Average fish MeHg Vapaa-ajan BAU Female 0.67722764
## 55 Average fish Omega3 Vapaa-ajan BAU Female 27.49754486
## 56 Average fish PFAS Vapaa-ajan BAU Female 27.57242657
## 57 Average fish TEQ Vapaa-ajan BAU Female 12.32895611
## 58 Average fish Vitamin D Vapaa-ajan BAU Female 0.41246317
## 59 Whitefish ALA Kasvatettu BAU Male 0.79278636
## 60 Whitefish DHA Kasvatettu BAU Male 1.40701723
## 61 Whitefish Fish Kasvatettu BAU Male 0.35711097
## 62 Whitefish MeHg Kasvatettu BAU Male 0.02363897
## 63 Whitefish Omega3 Kasvatettu BAU Male 3.57110972
## 64 Whitefish PFAS Kasvatettu BAU Male 3.14012142
## 65 Whitefish TEQ Kasvatettu BAU Male 2.05327336
## 66 Whitefish Vitamin D Kasvatettu BAU Male 0.05142398
## 67 Average fish ALA Kaupallinen BAU Male 1.31416838
## 68 Average fish DHA Kaupallinen BAU Male 4.83766330
## 69 Average fish Fish Kaupallinen BAU Male 1.90459185
## 70 Average fish MeHg Kaupallinen BAU Male 0.32835280
## 71 Average fish Omega3 Kaupallinen BAU Male 13.33214296
## 72 Average fish PFAS Kaupallinen BAU Male 13.36844925
## 73 Average fish TEQ Kaupallinen BAU Male 5.97767569
## 74 Average fish Vitamin D Kaupallinen BAU Male 0.19998214
## 75 Rainbow trout ALA Kirjolohi BAU Male 19.18102550
## 76 Rainbow trout DHA Kirjolohi BAU Male 30.18718567
## 77 Rainbow trout Fish Kirjolohi BAU Male 3.98773919
## 78 Rainbow trout Omega3 Kirjolohi BAU Male 71.77930542
## 79 Rainbow trout PFAS Kirjolohi BAU Male 35.06468916
## 80 Rainbow trout Vitamin D Kirjolohi BAU Male 0.20337470
## 81 Average fish ALA Muu tuonti BAU Male 11.49897331
## 82 Average fish DHA Muu tuonti BAU Male 42.32955391
## 83 Average fish Fish Muu tuonti BAU Male 16.66517870
## 84 Average fish MeHg Muu tuonti BAU Male 2.87308697
## 85 Average fish Omega3 Muu tuonti BAU Male 116.65625093
## 86 Average fish PFAS Muu tuonti BAU Male 116.97393090
## 87 Average fish TEQ Muu tuonti BAU Male 52.30466229
## 88 Average fish Vitamin D Muu tuonti BAU Male 1.74984376
## 89 Herring ALA Silakka BAU Male 1.65699491
## 90 Herring DHA Silakka BAU Male 5.58045413
## 91 Herring Fish Silakka BAU Male 0.95229593
## 92 Herring MeHg Silakka BAU Male 0.02553100
## 93 Herring Omega3 Silakka BAU Male 22.85510222
## 94 Herring PFAS Silakka BAU Male 1.57604976
## 95 Herring Vitamin D Silakka BAU Male 0.14855816
## 96 Rainbow trout ALA Tuontikirjolohi BAU Male 13.74163021
## 97 Rainbow trout DHA Tuontikirjolohi BAU Male 21.62664048
## 98 Rainbow trout Fish Tuontikirjolohi BAU Male 2.85688778
## 99 Rainbow trout Omega3 Tuontikirjolohi BAU Male 51.42398000
## 100 Rainbow trout PFAS Tuontikirjolohi BAU Male 25.12097134
## 101 Rainbow trout Vitamin D Tuontikirjolohi BAU Male 0.14570128
## 102 Salmon ALA Tuontilohi BAU Male 111.80906467
## 103 Salmon DHA Tuontilohi BAU Male 93.97018123
## 104 Salmon Fish Tuontilohi BAU Male 14.04636491
## 105 Salmon Omega3 Tuontilohi BAU Male 323.06639288
## 106 Salmon PFAS Tuontilohi BAU Male 123.51144243
## 107 Salmon TEQ Tuontilohi BAU Male 281.03662052
## 108 Salmon Vitamin D Tuontilohi BAU Male 0.94110645
## 109 Average fish ALA Vapaa-ajan BAU Male 4.06570842
## 110 Average fish DHA Vapaa-ajan BAU Male 14.96652085
## 111 Average fish Fish Vapaa-ajan BAU Male 5.89233104
## 112 Average fish MeHg Vapaa-ajan BAU Male 1.01584146
## 113 Average fish Omega3 Vapaa-ajan BAU Male 41.24631729
## 114 Average fish PFAS Vapaa-ajan BAU Male 41.35863985
## 115 Average fish TEQ Vapaa-ajan BAU Male 18.49343417
## 116 Average fish Vitamin D Vapaa-ajan BAU Male 0.61869476
#oprint(summary(exposure,"mean"))
#oprint(summary(fish_proportion,"mean"))
oprint(summary(incidence,"mean"))
## Response Age Adjust mean
## 1 Immunosuppression 0 - 4 BAU 4.8100
## 2 Loss in child's IQ points 0 - 4 BAU 1.1920
## 3 Sperm concentration 0 - 4 BAU 0.0140
## 4 Yes or no dental defect 0 - 4 BAU 0.0448
## 5 Immunosuppression 5 - 9 BAU 3.9500
## 6 Dioxin recommendation tolerable daily intake Undefined BAU 0.1100
## 7 Dioxin recommendation tolerable daily intake 2018 Undefined BAU 0.3200
## 8 PFAS TWI Undefined BAU 1.0000
## 9 Vitamin D recommendation Undefined BAU 0.2200
oprint(summary(PAF,"mean"))
## Exposure_agent Response Scenario
## 1 DHA Loss in child's IQ points BAU
## 2 MeHg Loss in child's IQ points BAU
## 3 TEQ Cancer morbidity BAU
## 4 TEQ Sperm concentration BAU
## 5 TEQ Yes or no dental defect BAU
## 6 PFAS Immunosuppression BAU
## 7 Fish All-cause mortality BAU
## 8 Omega3 Breast cancer BAU
## 9 Fish Depression BAU
## 10 Omega3 CHD2 mortality BAU
## 11 DHA Loss in child's IQ points BAU
## 12 MeHg Loss in child's IQ points BAU
## 13 TEQ Cancer morbidity BAU
## 14 TEQ Sperm concentration BAU
## 15 TEQ Yes or no dental defect BAU
## 16 PFAS Immunosuppression BAU
## 17 Fish All-cause mortality BAU
## 18 Omega3 Breast cancer BAU
## 19 Fish Depression BAU
## 20 Omega3 CHD2 mortality BAU
## 21 TEQ Cancer morbidity BAU
## 22 Fish All-cause mortality BAU
## 23 Omega3 Breast cancer BAU
## 24 Fish Depression BAU
## 25 Omega3 CHD2 mortality BAU
## 26 TEQ Cancer morbidity BAU
## 27 Fish All-cause mortality BAU
## 28 Omega3 Breast cancer BAU
## 29 Fish Depression BAU
## 30 Omega3 CHD2 mortality BAU
## 31 TEQ Cancer morbidity BAU
## 32 Fish All-cause mortality BAU
## 33 Omega3 Breast cancer BAU
## 34 Fish Depression BAU
## 35 Omega3 CHD2 mortality BAU
## 36 TEQ Cancer morbidity BAU
## 37 Fish All-cause mortality BAU
## 38 Omega3 Breast cancer BAU
## 39 Fish Depression BAU
## 40 Omega3 CHD2 mortality BAU
## 41 TEQ Cancer morbidity BAU
## 42 Fish All-cause mortality BAU
## 43 Omega3 Breast cancer BAU
## 44 Fish Depression BAU
## 45 Omega3 CHD2 mortality BAU
## 46 TEQ Cancer morbidity BAU
## 47 Fish All-cause mortality BAU
## 48 Omega3 Breast cancer BAU
## 49 Fish Depression BAU
## 50 Omega3 CHD2 mortality BAU
## 51 TEQ Cancer morbidity BAU
## 52 Fish All-cause mortality BAU
## 53 Omega3 Breast cancer BAU
## 54 Fish Depression BAU
## 55 Omega3 CHD2 mortality BAU
## 56 TEQ Cancer morbidity BAU
## 57 Fish All-cause mortality BAU
## 58 Omega3 Breast cancer BAU
## 59 Fish Depression BAU
## 60 Omega3 CHD2 mortality BAU
## 61 TEQ Cancer morbidity BAU
## 62 Fish All-cause mortality BAU
## 63 Omega3 Breast cancer BAU
## 64 Fish Depression BAU
## 65 Omega3 CHD2 mortality BAU
## 66 TEQ Cancer morbidity BAU
## 67 Fish All-cause mortality BAU
## 68 Omega3 Breast cancer BAU
## 69 Fish Depression BAU
## 70 Omega3 CHD2 mortality BAU
## 71 TEQ Cancer morbidity BAU
## 72 Fish All-cause mortality BAU
## 73 Omega3 Breast cancer BAU
## 74 Fish Depression BAU
## 75 Omega3 CHD2 mortality BAU
## 76 TEQ Cancer morbidity BAU
## 77 Fish All-cause mortality BAU
## 78 Omega3 Breast cancer BAU
## 79 Fish Depression BAU
## 80 Omega3 CHD2 mortality BAU
## 81 TEQ Cancer morbidity BAU
## 82 Fish All-cause mortality BAU
## 83 Omega3 Breast cancer BAU
## 84 Fish Depression BAU
## 85 Omega3 CHD2 mortality BAU
## 86 TEQ Cancer morbidity BAU
## 87 Fish All-cause mortality BAU
## 88 Omega3 Breast cancer BAU
## 89 Fish Depression BAU
## 90 Omega3 CHD2 mortality BAU
## 91 TEQ Cancer morbidity BAU
## 92 Fish All-cause mortality BAU
## 93 Omega3 Breast cancer BAU
## 94 Fish Depression BAU
## 95 Omega3 CHD2 mortality BAU
## 96 TEQ Cancer morbidity BAU
## 97 Fish All-cause mortality BAU
## 98 Omega3 Breast cancer BAU
## 99 Fish Depression BAU
## 100 Omega3 CHD2 mortality BAU
## 101 TEQ Cancer morbidity BAU
## 102 PFAS Immunosuppression BAU
## 103 Fish All-cause mortality BAU
## 104 Omega3 Breast cancer BAU
## 105 Fish Depression BAU
## 106 Omega3 CHD2 mortality BAU
## 107 TEQ Cancer morbidity BAU
## 108 PFAS Immunosuppression BAU
## 109 Fish All-cause mortality BAU
## 110 Omega3 Breast cancer BAU
## 111 Fish Depression BAU
## 112 Omega3 CHD2 mortality BAU
## 113 TEQ Cancer morbidity BAU
## 114 Fish All-cause mortality BAU
## 115 Omega3 Breast cancer BAU
## 116 Fish Depression BAU
## 117 Omega3 CHD2 mortality BAU
## 118 TEQ Cancer morbidity BAU
## 119 Fish All-cause mortality BAU
## 120 Omega3 Breast cancer BAU
## 121 Fish Depression BAU
## 122 Omega3 CHD2 mortality BAU
## 123 TEQ Cancer morbidity BAU
## 124 Fish All-cause mortality BAU
## 125 Omega3 Breast cancer BAU
## 126 Fish Depression BAU
## 127 Omega3 CHD2 mortality BAU
## 128 TEQ Cancer morbidity BAU
## 129 Fish All-cause mortality BAU
## 130 Omega3 Breast cancer BAU
## 131 Fish Depression BAU
## 132 Omega3 CHD2 mortality BAU
## 133 TEQ Cancer morbidity BAU
## 134 Fish All-cause mortality BAU
## 135 Omega3 Breast cancer BAU
## 136 Fish Depression BAU
## 137 Omega3 CHD2 mortality BAU
## 138 TEQ Cancer morbidity BAU
## 139 Fish All-cause mortality BAU
## 140 Omega3 Breast cancer BAU
## 141 Fish Depression BAU
## 142 Omega3 CHD2 mortality BAU
## 143 TEQ Cancer morbidity BAU
## 144 Fish All-cause mortality BAU
## 145 Omega3 Breast cancer BAU
## 146 Fish Depression BAU
## 147 Omega3 CHD2 mortality BAU
## 148 TEQ Cancer morbidity BAU
## 149 Fish All-cause mortality BAU
## 150 Omega3 Breast cancer BAU
## 151 Fish Depression BAU
## 152 Omega3 CHD2 mortality BAU
## 153 TEQ Cancer morbidity BAU
## 154 Fish All-cause mortality BAU
## 155 Omega3 Breast cancer BAU
## 156 Fish Depression BAU
## 157 Omega3 CHD2 mortality BAU
## 158 TEQ Cancer morbidity BAU
## 159 Fish All-cause mortality BAU
## 160 Omega3 Breast cancer BAU
## 161 Fish Depression BAU
## 162 Omega3 CHD2 mortality BAU
## 163 TEQ Cancer morbidity BAU
## 164 Fish All-cause mortality BAU
## 165 Omega3 Breast cancer BAU
## 166 Fish Depression BAU
## 167 Omega3 CHD2 mortality BAU
## 168 TEQ Cancer morbidity BAU
## 169 Fish All-cause mortality BAU
## 170 Omega3 Breast cancer BAU
## 171 Fish Depression BAU
## 172 Omega3 CHD2 mortality BAU
## 173 TEQ Cancer morbidity BAU
## 174 Fish All-cause mortality BAU
## 175 Omega3 Breast cancer BAU
## 176 Fish Depression BAU
## 177 Omega3 CHD2 mortality BAU
## 178 TEQ Cancer morbidity BAU
## 179 Fish All-cause mortality BAU
## 180 Omega3 Breast cancer BAU
## 181 Fish Depression BAU
## 182 Omega3 CHD2 mortality BAU
## 183 TEQ Cancer morbidity BAU
## 184 Fish All-cause mortality BAU
## 185 Omega3 Breast cancer BAU
## 186 Fish Depression BAU
## 187 Omega3 CHD2 mortality BAU
## 188 TEQ Cancer morbidity BAU
## 189 Fish All-cause mortality BAU
## 190 Omega3 Breast cancer BAU
## 191 Fish Depression BAU
## 192 Omega3 CHD2 mortality BAU
## 193 TEQ Cancer morbidity BAU
## 194 TEQ Dioxin recommendation tolerable daily intake BAU
## 195 TEQ Dioxin recommendation tolerable daily intake 2018 BAU
## 196 Vitamin D Vitamin D recommendation BAU
## 197 PFAS PFAS TWI BAU
## 198 Fish All-cause mortality BAU
## 199 Omega3 Breast cancer BAU
## 200 Fish Depression BAU
## 201 Omega3 CHD2 mortality BAU
## 202 TEQ Cancer morbidity BAU
## 203 TEQ Dioxin recommendation tolerable daily intake BAU
## 204 TEQ Dioxin recommendation tolerable daily intake 2018 BAU
## 205 Vitamin D Vitamin D recommendation BAU
## 206 PFAS PFAS TWI BAU
## 207 Fish All-cause mortality BAU
## 208 Omega3 Breast cancer BAU
## 209 Fish Depression BAU
## 210 Omega3 CHD2 mortality BAU
## Gender Age Adjust mean
## 1 Female 0 - 4 BAU -0.15625100
## 2 Female 0 - 4 BAU 0.22245956
## 3 Female 0 - 4 BAU 0.20518221
## 4 Female 0 - 4 BAU 0.80882827
## 5 Female 0 - 4 BAU 5.85966964
## 6 Female 0 - 4 BAU 0.01813625
## 7 Female 0 - 4 BAU -0.06412974
## 8 Female 0 - 4 BAU -0.19763757
## 9 Female 0 - 4 BAU -0.15262578
## 10 Female 0 - 4 BAU -0.15322440
## 11 Male 0 - 4 BAU -0.23437649
## 12 Male 0 - 4 BAU 0.33368933
## 13 Male 0 - 4 BAU 0.30777331
## 14 Male 0 - 4 BAU 1.21324240
## 15 Male 0 - 4 BAU 8.78950446
## 16 Male 0 - 4 BAU 0.02720438
## 17 Male 0 - 4 BAU -0.09463548
## 18 Male 0 - 4 BAU -0.28128638
## 19 Male 0 - 4 BAU -0.21996719
## 20 Male 0 - 4 BAU -0.15843589
## 21 Female 10 - 14 BAU 0.20518221
## 22 Female 10 - 14 BAU -0.06412974
## 23 Female 10 - 14 BAU -0.19763757
## 24 Female 10 - 14 BAU -0.15262578
## 25 Female 10 - 14 BAU -0.15322440
## 26 Male 10 - 14 BAU 0.30777331
## 27 Male 10 - 14 BAU -0.09463548
## 28 Male 10 - 14 BAU -0.28128638
## 29 Male 10 - 14 BAU -0.21996719
## 30 Male 10 - 14 BAU -0.15843589
## 31 Female 15 - 19 BAU 0.20518221
## 32 Female 15 - 19 BAU -0.06412974
## 33 Female 15 - 19 BAU -0.19763757
## 34 Female 15 - 19 BAU -0.15262578
## 35 Female 15 - 19 BAU -0.15322440
## 36 Male 15 - 19 BAU 0.30777331
## 37 Male 15 - 19 BAU -0.09463548
## 38 Male 15 - 19 BAU -0.28128638
## 39 Male 15 - 19 BAU -0.21996719
## 40 Male 15 - 19 BAU -0.15843589
## 41 Female 20 - 24 BAU 0.20518221
## 42 Female 20 - 24 BAU -0.06412974
## 43 Female 20 - 24 BAU -0.19763757
## 44 Female 20 - 24 BAU -0.15262578
## 45 Female 20 - 24 BAU -0.15322440
## 46 Male 20 - 24 BAU 0.30777331
## 47 Male 20 - 24 BAU -0.09463548
## 48 Male 20 - 24 BAU -0.28128638
## 49 Male 20 - 24 BAU -0.21996719
## 50 Male 20 - 24 BAU -0.15843589
## 51 Female 25 - 29 BAU 0.20518221
## 52 Female 25 - 29 BAU -0.06412974
## 53 Female 25 - 29 BAU -0.19763757
## 54 Female 25 - 29 BAU -0.15262578
## 55 Female 25 - 29 BAU -0.15322440
## 56 Male 25 - 29 BAU 0.30777331
## 57 Male 25 - 29 BAU -0.09463548
## 58 Male 25 - 29 BAU -0.28128638
## 59 Male 25 - 29 BAU -0.21996719
## 60 Male 25 - 29 BAU -0.15843589
## 61 Female 30 - 34 BAU 0.20518221
## 62 Female 30 - 34 BAU -0.06412974
## 63 Female 30 - 34 BAU -0.19763757
## 64 Female 30 - 34 BAU -0.15262578
## 65 Female 30 - 34 BAU -0.15322440
## 66 Male 30 - 34 BAU 0.30777331
## 67 Male 30 - 34 BAU -0.09463548
## 68 Male 30 - 34 BAU -0.28128638
## 69 Male 30 - 34 BAU -0.21996719
## 70 Male 30 - 34 BAU -0.15843589
## 71 Female 35 - 39 BAU 0.20518221
## 72 Female 35 - 39 BAU -0.06412974
## 73 Female 35 - 39 BAU -0.19763757
## 74 Female 35 - 39 BAU -0.15262578
## 75 Female 35 - 39 BAU -0.15322440
## 76 Male 35 - 39 BAU 0.30777331
## 77 Male 35 - 39 BAU -0.09463548
## 78 Male 35 - 39 BAU -0.28128638
## 79 Male 35 - 39 BAU -0.21996719
## 80 Male 35 - 39 BAU -0.15843589
## 81 Female 40 - 44 BAU 0.20518221
## 82 Female 40 - 44 BAU -0.06412974
## 83 Female 40 - 44 BAU -0.19763757
## 84 Female 40 - 44 BAU -0.15262578
## 85 Female 40 - 44 BAU -0.15322440
## 86 Male 40 - 44 BAU 0.30777331
## 87 Male 40 - 44 BAU -0.09463548
## 88 Male 40 - 44 BAU -0.28128638
## 89 Male 40 - 44 BAU -0.21996719
## 90 Male 40 - 44 BAU -0.15843589
## 91 Female 45 - 49 BAU 0.20518221
## 92 Female 45 - 49 BAU -0.06412974
## 93 Female 45 - 49 BAU -0.19763757
## 94 Female 45 - 49 BAU -0.15262578
## 95 Female 45 - 49 BAU -0.15322440
## 96 Male 45 - 49 BAU 0.30777331
## 97 Male 45 - 49 BAU -0.09463548
## 98 Male 45 - 49 BAU -0.28128638
## 99 Male 45 - 49 BAU -0.21996719
## 100 Male 45 - 49 BAU -0.15843589
## 101 Female 5 - 9 BAU 0.20518221
## 102 Female 5 - 9 BAU 0.02208490
## 103 Female 5 - 9 BAU -0.06412974
## 104 Female 5 - 9 BAU -0.19763757
## 105 Female 5 - 9 BAU -0.15262578
## 106 Female 5 - 9 BAU -0.15322440
## 107 Male 5 - 9 BAU 0.30777331
## 108 Male 5 - 9 BAU 0.03312735
## 109 Male 5 - 9 BAU -0.09463548
## 110 Male 5 - 9 BAU -0.28128638
## 111 Male 5 - 9 BAU -0.21996719
## 112 Male 5 - 9 BAU -0.15843589
## 113 Female 50 - 54 BAU 0.20518221
## 114 Female 50 - 54 BAU -0.06412974
## 115 Female 50 - 54 BAU -0.19763757
## 116 Female 50 - 54 BAU -0.15262578
## 117 Female 50 - 54 BAU -0.15322440
## 118 Male 50 - 54 BAU 0.30777331
## 119 Male 50 - 54 BAU -0.09463548
## 120 Male 50 - 54 BAU -0.28128638
## 121 Male 50 - 54 BAU -0.21996719
## 122 Male 50 - 54 BAU -0.15843589
## 123 Female 55 - 59 BAU 0.20518221
## 124 Female 55 - 59 BAU -0.06412974
## 125 Female 55 - 59 BAU -0.19763757
## 126 Female 55 - 59 BAU -0.15262578
## 127 Female 55 - 59 BAU -0.15322440
## 128 Male 55 - 59 BAU 0.30777331
## 129 Male 55 - 59 BAU -0.09463548
## 130 Male 55 - 59 BAU -0.28128638
## 131 Male 55 - 59 BAU -0.21996719
## 132 Male 55 - 59 BAU -0.15843589
## 133 Female 60 - 64 BAU 0.20518221
## 134 Female 60 - 64 BAU -0.06412974
## 135 Female 60 - 64 BAU -0.19763757
## 136 Female 60 - 64 BAU -0.15262578
## 137 Female 60 - 64 BAU -0.15322440
## 138 Male 60 - 64 BAU 0.30777331
## 139 Male 60 - 64 BAU -0.09463548
## 140 Male 60 - 64 BAU -0.28128638
## 141 Male 60 - 64 BAU -0.21996719
## 142 Male 60 - 64 BAU -0.15843589
## 143 Female 65 - 69 BAU 0.20518221
## 144 Female 65 - 69 BAU -0.06412974
## 145 Female 65 - 69 BAU -0.19763757
## 146 Female 65 - 69 BAU -0.15262578
## 147 Female 65 - 69 BAU -0.15322440
## 148 Male 65 - 69 BAU 0.30777331
## 149 Male 65 - 69 BAU -0.09463548
## 150 Male 65 - 69 BAU -0.28128638
## 151 Male 65 - 69 BAU -0.21996719
## 152 Male 65 - 69 BAU -0.15843589
## 153 Female 70 - 74 BAU 0.20518221
## 154 Female 70 - 74 BAU -0.06412974
## 155 Female 70 - 74 BAU -0.19763757
## 156 Female 70 - 74 BAU -0.15262578
## 157 Female 70 - 74 BAU -0.15322440
## 158 Male 70 - 74 BAU 0.30777331
## 159 Male 70 - 74 BAU -0.09463548
## 160 Male 70 - 74 BAU -0.28128638
## 161 Male 70 - 74 BAU -0.21996719
## 162 Male 70 - 74 BAU -0.15843589
## 163 Female 75 - 79 BAU 0.20518221
## 164 Female 75 - 79 BAU -0.06412974
## 165 Female 75 - 79 BAU -0.19763757
## 166 Female 75 - 79 BAU -0.15262578
## 167 Female 75 - 79 BAU -0.15322440
## 168 Male 75 - 79 BAU 0.30777331
## 169 Male 75 - 79 BAU -0.09463548
## 170 Male 75 - 79 BAU -0.28128638
## 171 Male 75 - 79 BAU -0.21996719
## 172 Male 75 - 79 BAU -0.15843589
## 173 Female 80 - 84 BAU 0.20518221
## 174 Female 80 - 84 BAU -0.06412974
## 175 Female 80 - 84 BAU -0.19763757
## 176 Female 80 - 84 BAU -0.15262578
## 177 Female 80 - 84 BAU -0.15322440
## 178 Male 80 - 84 BAU 0.30777331
## 179 Male 80 - 84 BAU -0.09463548
## 180 Male 80 - 84 BAU -0.28128638
## 181 Male 80 - 84 BAU -0.21996719
## 182 Male 80 - 84 BAU -0.15843589
## 183 Female 85+ BAU 0.20518221
## 184 Female 85+ BAU -0.06412974
## 185 Female 85+ BAU -0.19763757
## 186 Female 85+ BAU -0.15262578
## 187 Female 85+ BAU -0.15322440
## 188 Male 85+ BAU 0.30777331
## 189 Male 85+ BAU -0.09463548
## 190 Male 85+ BAU -0.28128638
## 191 Male 85+ BAU -0.21996719
## 192 Male 85+ BAU -0.15843589
## 193 Female Undefined BAU 0.20518221
## 194 Female Undefined BAU 1.00000000
## 195 Female Undefined BAU 1.00000000
## 196 Female Undefined BAU 1.00000000
## 197 Female Undefined BAU 0.00000000
## 198 Female Undefined BAU -0.06412974
## 199 Female Undefined BAU -0.19763757
## 200 Female Undefined BAU -0.15262578
## 201 Female Undefined BAU -0.15322440
## 202 Male Undefined BAU 0.30777331
## 203 Male Undefined BAU 1.00000000
## 204 Male Undefined BAU 1.00000000
## 205 Male Undefined BAU 1.00000000
## 206 Male Undefined BAU 1.00000000
## 207 Male Undefined BAU -0.09463548
## 208 Male Undefined BAU -0.28128638
## 209 Male Undefined BAU -0.21996719
## 210 Male Undefined BAU -0.15843589
oprint(summary(population,"mean"))
## Gender Age mean
## 1 Female 0 - 4 125040
## 2 Male 0 - 4 130884
## 3 Female 5 - 9 149633
## 4 Male 5 - 9 156654
## 5 Female 10 - 14 151113
## 6 Male 10 - 14 157712
## 7 Female 15 - 19 144441
## 8 Male 15 - 19 152230
## 9 Female 20 - 24 152265
## 10 Male 20 - 24 161679
## 11 Female 25 - 29 172593
## 12 Male 25 - 29 183092
## 13 Female 30 - 34 169653
## 14 Male 30 - 34 181115
## 15 Female 35 - 39 174660
## 16 Male 35 - 39 186122
## 17 Female 40 - 44 168547
## 18 Male 40 - 44 177928
## 19 Female 45 - 49 154391
## 20 Male 45 - 49 159982
## 21 Female 50 - 54 176612
## 22 Male 50 - 54 179182
## 23 Female 55 - 59 185152
## 24 Male 55 - 59 183719
## 25 Female 60 - 64 183336
## 26 Male 60 - 64 176283
## 27 Female 65 - 69 185685
## 28 Male 65 - 69 171275
## 29 Female 70 - 74 186034
## 30 Male 70 - 74 163697
## 31 Female 75 - 79 118190
## 32 Male 75 - 79 93987
## 33 Female 80 - 84 96256
## 34 Male 80 - 84 65140
## 35 Female 85+ 103429
## 36 Male 85+ 47581
## 37 Female Undefined 2797030
## 38 Male Undefined 2728262
oprint(summary(RR,"mean"))
## Exposure_agent Response ER_function Scaling Scenario Gender
## 1 Fish All-cause mortality RR None BAU Female
## 2 Omega3 Breast cancer RR None BAU Female
## 3 Fish Depression RR None BAU Female
## 4 Omega3 CHD2 mortality Relative Hill None BAU Female
## 5 Fish All-cause mortality RR None BAU Male
## 6 Omega3 Breast cancer RR None BAU Male
## 7 Fish Depression RR None BAU Male
## 8 Omega3 CHD2 mortality Relative Hill None BAU Male
## 9 Fish All-cause mortality RR None BAU Female
## 10 Omega3 Breast cancer RR None BAU Female
## 11 Fish Depression RR None BAU Female
## 12 Omega3 CHD2 mortality Relative Hill None BAU Female
## 13 Fish All-cause mortality RR None BAU Male
## 14 Omega3 Breast cancer RR None BAU Male
## 15 Fish Depression RR None BAU Male
## 16 Omega3 CHD2 mortality Relative Hill None BAU Male
## 17 Fish All-cause mortality RR None BAU Female
## 18 Omega3 Breast cancer RR None BAU Female
## 19 Fish Depression RR None BAU Female
## 20 Omega3 CHD2 mortality Relative Hill None BAU Female
## 21 Fish All-cause mortality RR None BAU Male
## 22 Omega3 Breast cancer RR None BAU Male
## 23 Fish Depression RR None BAU Male
## 24 Omega3 CHD2 mortality Relative Hill None BAU Male
## 25 Fish All-cause mortality RR None BAU Female
## 26 Omega3 Breast cancer RR None BAU Female
## 27 Fish Depression RR None BAU Female
## 28 Omega3 CHD2 mortality Relative Hill None BAU Female
## 29 Fish All-cause mortality RR None BAU Male
## 30 Omega3 Breast cancer RR None BAU Male
## 31 Fish Depression RR None BAU Male
## 32 Omega3 CHD2 mortality Relative Hill None BAU Male
## 33 Fish All-cause mortality RR None BAU Female
## 34 Omega3 Breast cancer RR None BAU Female
## 35 Fish Depression RR None BAU Female
## 36 Omega3 CHD2 mortality Relative Hill None BAU Female
## 37 Fish All-cause mortality RR None BAU Male
## 38 Omega3 Breast cancer RR None BAU Male
## 39 Fish Depression RR None BAU Male
## 40 Omega3 CHD2 mortality Relative Hill None BAU Male
## 41 Fish All-cause mortality RR None BAU Female
## 42 Omega3 Breast cancer RR None BAU Female
## 43 Fish Depression RR None BAU Female
## 44 Omega3 CHD2 mortality Relative Hill None BAU Female
## 45 Fish All-cause mortality RR None BAU Male
## 46 Omega3 Breast cancer RR None BAU Male
## 47 Fish Depression RR None BAU Male
## 48 Omega3 CHD2 mortality Relative Hill None BAU Male
## 49 Fish All-cause mortality RR None BAU Female
## 50 Omega3 Breast cancer RR None BAU Female
## 51 Fish Depression RR None BAU Female
## 52 Omega3 CHD2 mortality Relative Hill None BAU Female
## 53 Fish All-cause mortality RR None BAU Male
## 54 Omega3 Breast cancer RR None BAU Male
## 55 Fish Depression RR None BAU Male
## 56 Omega3 CHD2 mortality Relative Hill None BAU Male
## 57 Fish All-cause mortality RR None BAU Female
## 58 Omega3 Breast cancer RR None BAU Female
## 59 Fish Depression RR None BAU Female
## 60 Omega3 CHD2 mortality Relative Hill None BAU Female
## 61 Fish All-cause mortality RR None BAU Male
## 62 Omega3 Breast cancer RR None BAU Male
## 63 Fish Depression RR None BAU Male
## 64 Omega3 CHD2 mortality Relative Hill None BAU Male
## 65 Fish All-cause mortality RR None BAU Female
## 66 Omega3 Breast cancer RR None BAU Female
## 67 Fish Depression RR None BAU Female
## 68 Omega3 CHD2 mortality Relative Hill None BAU Female
## 69 Fish All-cause mortality RR None BAU Male
## 70 Omega3 Breast cancer RR None BAU Male
## 71 Fish Depression RR None BAU Male
## 72 Omega3 CHD2 mortality Relative Hill None BAU Male
## 73 Fish All-cause mortality RR None BAU Female
## 74 Omega3 Breast cancer RR None BAU Female
## 75 Fish Depression RR None BAU Female
## 76 Omega3 CHD2 mortality Relative Hill None BAU Female
## 77 Fish All-cause mortality RR None BAU Male
## 78 Omega3 Breast cancer RR None BAU Male
## 79 Fish Depression RR None BAU Male
## 80 Omega3 CHD2 mortality Relative Hill None BAU Male
## 81 Fish All-cause mortality RR None BAU Female
## 82 Omega3 Breast cancer RR None BAU Female
## 83 Fish Depression RR None BAU Female
## 84 Omega3 CHD2 mortality Relative Hill None BAU Female
## 85 Fish All-cause mortality RR None BAU Male
## 86 Omega3 Breast cancer RR None BAU Male
## 87 Fish Depression RR None BAU Male
## 88 Omega3 CHD2 mortality Relative Hill None BAU Male
## 89 Fish All-cause mortality RR None BAU Female
## 90 Omega3 Breast cancer RR None BAU Female
## 91 Fish Depression RR None BAU Female
## 92 Omega3 CHD2 mortality Relative Hill None BAU Female
## 93 Fish All-cause mortality RR None BAU Male
## 94 Omega3 Breast cancer RR None BAU Male
## 95 Fish Depression RR None BAU Male
## 96 Omega3 CHD2 mortality Relative Hill None BAU Male
## 97 Fish All-cause mortality RR None BAU Female
## 98 Omega3 Breast cancer RR None BAU Female
## 99 Fish Depression RR None BAU Female
## 100 Omega3 CHD2 mortality Relative Hill None BAU Female
## 101 Fish All-cause mortality RR None BAU Male
## 102 Omega3 Breast cancer RR None BAU Male
## 103 Fish Depression RR None BAU Male
## 104 Omega3 CHD2 mortality Relative Hill None BAU Male
## 105 Fish All-cause mortality RR None BAU Female
## 106 Omega3 Breast cancer RR None BAU Female
## 107 Fish Depression RR None BAU Female
## 108 Omega3 CHD2 mortality Relative Hill None BAU Female
## 109 Fish All-cause mortality RR None BAU Male
## 110 Omega3 Breast cancer RR None BAU Male
## 111 Fish Depression RR None BAU Male
## 112 Omega3 CHD2 mortality Relative Hill None BAU Male
## 113 Fish All-cause mortality RR None BAU Female
## 114 Omega3 Breast cancer RR None BAU Female
## 115 Fish Depression RR None BAU Female
## 116 Omega3 CHD2 mortality Relative Hill None BAU Female
## 117 Fish All-cause mortality RR None BAU Male
## 118 Omega3 Breast cancer RR None BAU Male
## 119 Fish Depression RR None BAU Male
## 120 Omega3 CHD2 mortality Relative Hill None BAU Male
## 121 Fish All-cause mortality RR None BAU Female
## 122 Omega3 Breast cancer RR None BAU Female
## 123 Fish Depression RR None BAU Female
## 124 Omega3 CHD2 mortality Relative Hill None BAU Female
## 125 Fish All-cause mortality RR None BAU Male
## 126 Omega3 Breast cancer RR None BAU Male
## 127 Fish Depression RR None BAU Male
## 128 Omega3 CHD2 mortality Relative Hill None BAU Male
## 129 Fish All-cause mortality RR None BAU Female
## 130 Omega3 Breast cancer RR None BAU Female
## 131 Fish Depression RR None BAU Female
## 132 Omega3 CHD2 mortality Relative Hill None BAU Female
## 133 Fish All-cause mortality RR None BAU Male
## 134 Omega3 Breast cancer RR None BAU Male
## 135 Fish Depression RR None BAU Male
## 136 Omega3 CHD2 mortality Relative Hill None BAU Male
## 137 Fish All-cause mortality RR None BAU Female
## 138 Omega3 Breast cancer RR None BAU Female
## 139 Fish Depression RR None BAU Female
## 140 Omega3 CHD2 mortality Relative Hill None BAU Female
## 141 Fish All-cause mortality RR None BAU Male
## 142 Omega3 Breast cancer RR None BAU Male
## 143 Fish Depression RR None BAU Male
## 144 Omega3 CHD2 mortality Relative Hill None BAU Male
## 145 Fish All-cause mortality RR None BAU Female
## 146 Omega3 Breast cancer RR None BAU Female
## 147 Fish Depression RR None BAU Female
## 148 Omega3 CHD2 mortality Relative Hill None BAU Female
## 149 Fish All-cause mortality RR None BAU Male
## 150 Omega3 Breast cancer RR None BAU Male
## 151 Fish Depression RR None BAU Male
## 152 Omega3 CHD2 mortality Relative Hill None BAU Male
## Age mean
## 1 0 - 4 0.9358703
## 2 0 - 4 0.8023624
## 3 0 - 4 0.8473742
## 4 0 - 4 0.8467756
## 5 0 - 4 0.9053645
## 6 0 - 4 0.7187136
## 7 0 - 4 0.7800328
## 8 0 - 4 0.8415641
## 9 10 - 14 0.9358703
## 10 10 - 14 0.8023624
## 11 10 - 14 0.8473742
## 12 10 - 14 0.8467756
## 13 10 - 14 0.9053645
## 14 10 - 14 0.7187136
## 15 10 - 14 0.7800328
## 16 10 - 14 0.8415641
## 17 15 - 19 0.9358703
## 18 15 - 19 0.8023624
## 19 15 - 19 0.8473742
## 20 15 - 19 0.8467756
## 21 15 - 19 0.9053645
## 22 15 - 19 0.7187136
## 23 15 - 19 0.7800328
## 24 15 - 19 0.8415641
## 25 20 - 24 0.9358703
## 26 20 - 24 0.8023624
## 27 20 - 24 0.8473742
## 28 20 - 24 0.8467756
## 29 20 - 24 0.9053645
## 30 20 - 24 0.7187136
## 31 20 - 24 0.7800328
## 32 20 - 24 0.8415641
## 33 25 - 29 0.9358703
## 34 25 - 29 0.8023624
## 35 25 - 29 0.8473742
## 36 25 - 29 0.8467756
## 37 25 - 29 0.9053645
## 38 25 - 29 0.7187136
## 39 25 - 29 0.7800328
## 40 25 - 29 0.8415641
## 41 30 - 34 0.9358703
## 42 30 - 34 0.8023624
## 43 30 - 34 0.8473742
## 44 30 - 34 0.8467756
## 45 30 - 34 0.9053645
## 46 30 - 34 0.7187136
## 47 30 - 34 0.7800328
## 48 30 - 34 0.8415641
## 49 35 - 39 0.9358703
## 50 35 - 39 0.8023624
## 51 35 - 39 0.8473742
## 52 35 - 39 0.8467756
## 53 35 - 39 0.9053645
## 54 35 - 39 0.7187136
## 55 35 - 39 0.7800328
## 56 35 - 39 0.8415641
## 57 40 - 44 0.9358703
## 58 40 - 44 0.8023624
## 59 40 - 44 0.8473742
## 60 40 - 44 0.8467756
## 61 40 - 44 0.9053645
## 62 40 - 44 0.7187136
## 63 40 - 44 0.7800328
## 64 40 - 44 0.8415641
## 65 45 - 49 0.9358703
## 66 45 - 49 0.8023624
## 67 45 - 49 0.8473742
## 68 45 - 49 0.8467756
## 69 45 - 49 0.9053645
## 70 45 - 49 0.7187136
## 71 45 - 49 0.7800328
## 72 45 - 49 0.8415641
## 73 5 - 9 0.9358703
## 74 5 - 9 0.8023624
## 75 5 - 9 0.8473742
## 76 5 - 9 0.8467756
## 77 5 - 9 0.9053645
## 78 5 - 9 0.7187136
## 79 5 - 9 0.7800328
## 80 5 - 9 0.8415641
## 81 50 - 54 0.9358703
## 82 50 - 54 0.8023624
## 83 50 - 54 0.8473742
## 84 50 - 54 0.8467756
## 85 50 - 54 0.9053645
## 86 50 - 54 0.7187136
## 87 50 - 54 0.7800328
## 88 50 - 54 0.8415641
## 89 55 - 59 0.9358703
## 90 55 - 59 0.8023624
## 91 55 - 59 0.8473742
## 92 55 - 59 0.8467756
## 93 55 - 59 0.9053645
## 94 55 - 59 0.7187136
## 95 55 - 59 0.7800328
## 96 55 - 59 0.8415641
## 97 60 - 64 0.9358703
## 98 60 - 64 0.8023624
## 99 60 - 64 0.8473742
## 100 60 - 64 0.8467756
## 101 60 - 64 0.9053645
## 102 60 - 64 0.7187136
## 103 60 - 64 0.7800328
## 104 60 - 64 0.8415641
## 105 65 - 69 0.9358703
## 106 65 - 69 0.8023624
## 107 65 - 69 0.8473742
## 108 65 - 69 0.8467756
## 109 65 - 69 0.9053645
## 110 65 - 69 0.7187136
## 111 65 - 69 0.7800328
## 112 65 - 69 0.8415641
## 113 70 - 74 0.9358703
## 114 70 - 74 0.8023624
## 115 70 - 74 0.8473742
## 116 70 - 74 0.8467756
## 117 70 - 74 0.9053645
## 118 70 - 74 0.7187136
## 119 70 - 74 0.7800328
## 120 70 - 74 0.8415641
## 121 75 - 79 0.9358703
## 122 75 - 79 0.8023624
## 123 75 - 79 0.8473742
## 124 75 - 79 0.8467756
## 125 75 - 79 0.9053645
## 126 75 - 79 0.7187136
## 127 75 - 79 0.7800328
## 128 75 - 79 0.8415641
## 129 80 - 84 0.9358703
## 130 80 - 84 0.8023624
## 131 80 - 84 0.8473742
## 132 80 - 84 0.8467756
## 133 80 - 84 0.9053645
## 134 80 - 84 0.7187136
## 135 80 - 84 0.7800328
## 136 80 - 84 0.8415641
## 137 85+ 0.9358703
## 138 85+ 0.8023624
## 139 85+ 0.8473742
## 140 85+ 0.8467756
## 141 85+ 0.9053645
## 142 85+ 0.7187136
## 143 85+ 0.7800328
## 144 85+ 0.8415641
## 145 Undefined 0.9358703
## 146 Undefined 0.8023624
## 147 Undefined 0.8473742
## 148 Undefined 0.8467756
## 149 Undefined 0.9053645
## 150 Undefined 0.7187136
## 151 Undefined 0.7800328
## 152 Undefined 0.8415641
###################
# Graphs
trim <- function(ova) return(oapply(ova, NULL, mean, "Iter")@output)
ggplot(amount@output, aes(x=Gender, weight=amountResult, fill=Kala))+geom_bar()+
labs(
title="Kalansyönti Suomessa ikäryhmittäin",
y="Syönti (g/d)"
)
ggsave("Kalansyönti Suomessa ikäryhmittäin.svg")
## Saving 7 x 5 in image
plot_ly(trim(total_amount), x=~Kala, y=~total_amountResult, color=~Kala, type="bar") %>%
layout(yaxis=list(title="Kalan kokonaiskulutus Suomessa (milj kg /a)"), barmode="stack")
## Warning: `arrange_()` is deprecated as of dplyr 0.7.0.
## Please use `arrange()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.
plot_ly(trim(conc_vit), x=~Nutrient, y=~conc_vitResult, color=~Kala, type="scatter", mode="markers") %>%
layout(yaxis=list(title="Concentrations of nutrients (mg or ug /g)"))
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
ggplot(conc_pfas@output, aes(x=Fish, y=conc_pfasResult))+geom_point()+ # Inputed data for missing species.
coord_flip()
ggplot(conc_pfas@output, aes(x=conc_pfasResult, color=Compound, linetype=Area))+stat_ecdf()+
scale_x_log10()+
stat_ecdf(data=conc_eukalat@output, aes(x=conc_eukalatResult))+
scale_linetype_manual(values=c("dotted","solid","twodash"))+
labs(
title="PFAS concentration in fishes in Finland",
x="PFAS concentration (ng/g fresh weight)",
y="Cumulative probability"
)
# The code may produce some negative values, which are removed from the graph
ggsave("PFAS-pitoisuus kalassa Suomessa.svg")
## Saving 7 x 5 in image
ggplot(conc@output, aes(x=concResult, colour=Fish))+stat_ecdf()+
facet_wrap(~Compound, scales="free_x")
#ggplot(oapply(expo_dir, cols=c("Iter"),FUN=mean)@output,
# aes(x=Gender, weight=expo_dirResult,fill=Fish))+geom_bar()+
# facet_wrap(~Compound, scales="free_y")+
# labs(title="Eri yhdisteiden saanti kalasta")
#ggsave("Yhdisteiden saanti kalasta Suomessa.svg")
plot_ly(trim(exposure), x=~Gender, y=~exposureResult, color=~Exposure_agent, text=~Exposure_agent, type="bar") %>%
layout(yaxis=list(title="Exposure to nutrients (g or ug /d)"))
ggplot(trim(exposure), aes(x=Gender, weight=exposureResult, fill=Gender))+geom_bar()+
facet_wrap(~Exposure_agent, scales="free_y")+
labs(
title="Exposure to compounds",
y="(omega: mg/d; vit D: ug/d, PFAS: ng/d)"
)
cat("Kalaperäisiä tautitaakkoja Suomessa\n")
## Kalaperäisiä tautitaakkoja Suomessa
if(openv$N>1) {
tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Response")))
tmp <- data.frame(
Altiste = tmp$Exposure_agent,
Vaikutus = tmp$Response,
Keskiarvo = as.character(signif(tmp$mean,2)),
"95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
Keskihajonta = signif(tmp$sd,2)
)#[rev(match(lev, tmp$Exposure_agent)),]
oprint(tmp)
tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Exposure_agent")))
tmp <- data.frame(
Terveysvaikutus = tmp$Response,
Keskiarvo = signif(tmp$mean,2),
"95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
Keskihajonta = signif(tmp$sd,2)
)
oprint(tmp)
}
ggplot(trim(BoDattr), aes(x=Exposure_agent, weight=BoDattrResult, fill=Response))+geom_bar()
ggplot(trim(BoDattr), aes(x=Response, weight=BoDattrResult, fill=Exposure_agent))+geom_bar()
ggplot(trim(BoDattr), aes(x=Age, weight=BoDattrResult, fill=Response))+geom_bar(position="stack")
plot_ly(trim(BoDattr), x=~Exposure_agent, y=~BoDattrResult, color=~Response, text=~paste(Age, Exposure_agent, sep=": "), type="bar") %>%
layout(yaxis=list(title="Disease burden (DALY /a); CHD2=coronary heart disease"), barmode="stack")
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
################ Insight network
gr <- scrape(type="assessment")
objects.latest("Op_en3861", "makeGraph") # [[Insight network]]
## Loading objects:
## makeGraph
gr <- makeGraph(gr)
## Loading required package: DiagrammeR
## Loading objects:
## formatted
## Loading objects:
## chooseGr
#export_graph(gr, "ruori.svg")
#render_graph(gr) # Does not work: Error in generate_dot(graph) : object 'attribute' not found
##################### Diagnostics
objects.latest("Op_en6007", code_name="diagnostics")
## Loading objects:
## showind
## binoptest
## showLoctable
## ovashapetest
showLoctable()
showind()
## Hg is not an ovariable.
## subgrouping is not an ovariable.
## sumExposcen is not an ovariable.
## mc2d is not an ovariable.
## mc2dparam is not an ovariable.
## TEFversion is not an ovariable.